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Abstract 

Using data acquired with the CLEO-c detector at the CESR e^e^ colHder, we measure 
branching fractions for J/V', ip{2S), and ■i/'(3770) decays to 'jn^, 7ry, and 777'. Defining Rn = 
B[ip{nS) j'r]]/B[ip{nS) jv'], we obtain Ri = (21.1 ± 0.9)% and, unexpectedly, an order of 
magnitude smaller limit, R2 < 1.8% at 90% C.L. We also use J/il^^j-q' events to determine 
branching fractions of improved precision for the five most copious 77' decay modes. 

PACS numbers: 13.20.Gd,13.20.Jf 
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In the conventional view of quantum chromodynamics (QCD), the mass eigenstates J/V", 
ip{2S), and '?/'(3770) represent the IS', 25", and ID cc states, respectively (with some evidence 
of modest 2S — ID mixing py]). Most hadronic decays of these states that are not transitions 
to lower-lying charmonia or decays to open charm can be described at tree level by cc 
annihilation into either three gluons (ggg) or a photon plus two gluons (jgg) . However, some 
final states can be reached by less common routes, and their study can lend experimental 
constraints on the relevant QCD predictions. In this article we describe measurements of 
branching fractions for each of J/V', '0(25'), and ?/'(3770) decays to three such final states: 
77r°, 7?7, and 717'. In the hb system, limits of ;B(T(15') 'Ji]^'^) < 1.0(1.9) x 10^^ have been 
set |2| and are smaller than expected. Comparable studies in the charmonium sector are 
warranted. 

Vector charmonium decays to a photon and a hght pseudoscalar (ip^nS) 7P) can be 
described by a variety of mechanisms at the parton level, which we label {i) — {vi). When P = 
1] or rj', the primary source is generally thought to be {i) tplnS) —>■ '~fgg, which then fragments 
to exclusive final states in a flavor-blind manner. The r] and rj' mesons are commonly 
understood as mixtures of the pure SU(3)-fiavor octet [{uu + dd—2ss)/^/Q] and singlet [{uu + 
dd + ss) /a/3] states and a small gluonium component [1, 0]. The flavor content of each r)^'^ 
mass eigenstate can be quantified with a mixing angle, the value of which becomes manifest 
in ratios of branching fractions for various radiative decays involving an rj or r]' js], 0, 0] ■ 
Under the assumption that all ip{nS) — >■ 777^'^ final states occur through (i), the r]—T] ' mixing 
angle can be extracted from a measurement of Rn = B[ip{nS) 'yri]/B[ip{nS) •yr]'] and 
compared with values obtained from other decays [1, [1] . 

Alternate cc annihilation mechanisms for 7P final states are mediated either by three 
gluons or a virtual photon: (m) ip{nS) ggg qqjisi (where 7fsr represents final state 
radiation off one of the quarks) or {in) ip^nS) 7* —» qq7fsr- It has also been suggested 
that an Ml transition to rjc followed by r]c-r]^'^ mixing 

7Vr* IV could 

contribute. The vector dominance model (VDM) can be used to predict [ll| (v) th(2S) — 
J — ^ 777. For all these processes Ri ~ R2 is expected. Previous measurements |l2|] yield 
Ri = (20.8 ± 2.4)%, consistent with the dominance of (i) with an 77-77' mixing angle in the 
expected range [l2|. However, rate determinations for ip{2S) decays have yielded R2 < 66% 



at 90% confidence level (C.L.) pjj and hence are not yet precise enough to confirm a value 
comparable to -Ri. No measurements for ip(377Q) — > 777^'^ have yet been reported. 

When P = 71^, processes (i) and {iv) violate isospin conservation, and are therefore 
suppressed relative to process (iii) or (v) above, or (vi) ggg 2{qq) — > p°*7r°, p°* — >• 7 via 
VDM. Ref. [ill (CZ) finds process (ii) to be negligible, but [iii) and (vi) to be of comparable 
magnitude and fully coherent. Updating the CZ VDM calculation to current experimental 
branching fractions and widths I4I gives B(J/ijj — > 77?'') 6 x 10~^, nearly double the PDG 
average measurement of (3.31q;4) x 10~^ [l^l- A similar disparity was measured in Ref. [3] 
for 7*7 — i> tt" for spacelike non-asymptotic momentum transfers in the range = 1.5 — 9 
GeV^, where the CZ prediction was found to significantly overshoot the data. Further 
experimental precision on B{J/ip 77?°) would allow a mo re p recise indirect measurement 
of the 7* — 77r° vertex of process (iii) for timelike photons [l5]. 

For this measurement, events were acquired at the CESR e~'"e~ collider with the CLEO 
detector [3], mostly in the CLEO-c configuration (5% of 'M2S) data were collected with 
CLEO III). The data samples comprise (27.4 ± 0.6) x 10^ [l7| produced '(p{2S) mesons, of 
which (9.589 ± 0.020 ± 0.070) x 10^ Q decay into TT+n-J/ip (our source for J/ip decays), 
and 814 pb'^ at = 3.773 GeV, corresponding to (5.3 ± 0.1) x 10^ [Ts) produced ?/;(3770) 
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mesons. We also utilize 20.7 pb~^ of data taken at y/s = 3.671 GeV, just below the 4'{2S) 
peak, for a "continuum" subtraction of backgrounds. 

We select events in the decay modes shown in Table [H modes that are not shown have sen- 
sitivity inferior to those we employ. Every particle in each mode's decay chain must be found. 
Each such particle is constrained to originate from a single point (vertex) consistent with the 
measured beam spot. We then constrain the sum of all four- momenta to the known 'ip{2S) 



mass 



12| and initial e~^e~ three-momentum. The vertex and full event four-momentum 
kinematic fits must satisfy quality requirements of Xy/d.o.f. < 10 and xi/d-o.f. < 10, 
respectively, which typically retain 90-99% of signal events. For ip{2S) —>■ 77r°, tighter re- 
strictions of Xy/d.o.f. < 3 and xl/d.o.f. < 3 suppress QED background. Further selections 
are based on four- momenta from the fit. 

Masses of final state or intermediate mesons must be found in the signal or sideband 
windows given in Table [Tll Note that different windows apply depending on the t]^''' decay 
mode because mass resolution depends upon the final state. The tracking system must find 
zero, two, or four tracks of net charge zero for the n^/rj/r]' decay products, and, for J/ijj 
modes, exactly two additional oppositely-charged particles for the ?/'(25')-to-J/'i/'-transition. 
Photon candidate showers must lie in the barrel region of the calorimeter (| cos^| < 0.75) 
where the energy resolution and finding-efficiency are most favorable and well-modeled. 
Photon candidates must also have energy > 37 MeV, lie more than 30 cm from the 
center of any shower associated with one of the charged pions, and not align with the initial 
three-momentum of any vr^ candidate within 100 mrad. A photon candidate is included 
in the decay chain only if all more energetic ones are also used. Photon pairs from 7r° 
or 1] are constrained to the appropriate mass [l2|, except in 77 ^ Sn^ and 77' 7r°7r°?7(77); 
these final states are only required to have six photons with a combined invariant mass in the 
corresponding mass window in Table HTl Photons in rj^''^ — > 77? "'"tt" and rj' — > TT~^TT~ri('yTr~^TT~) 
must have > 100 MeV to suppress feedacross from other rj and 1]' decays; photons in 
ipi^S) 77r° and ^/'(3770) •yn^ must have E^^ > 300 MeV to suppress e~^e^ 37. 

For ^/'(3770) modes, where signal yields are expected to be small, backgrounds larger, and 
systematic considerations much less important, we take measures to enhance efficiency and 
suppress backgrounds at the expense of systematic precision. The photon-finding criteria 
are loosened so as to increase efficiency: barrel and endcap showers are used (| cos^|<0.83 
or 0.85<| cos^^|<0.93), and showers between 15 and 30 cm from a shower associated with a 
charged track are accepted if they have a photon-like lateral shape. In addition, the four- 
momentum constrained fit quality is tightened to x|;/d.o.f. < 3 to suppress backgrounds. 

Efficiencies for signal and feedacross from other modes considered here are modeled with 
Monte Carlo (MC) samples that were generated using the EvtGen event generator M, fed 
through a GEANT-based jloj detector simulation, and subjected to the same event selection 
criteria as the data. We model 77 — 77 and 37r° as phase-space-like decays. The decay 
T] 77r+7r~ is simulated as mediated by a p° — ^ tt+tt^, weighted with a factor E^, where E^ 
is the photon energy in the t] center-of-mass system, which results in the excellent agreement 



with the measured dipion mass distribution seen in Ref. [21[. The decay r]' 77? "•"vr is 
handled similarly, but using a factor of E^, and MC-data agreement is shown in Ref. [22 . 
We generate rj — ^ 7r+7r~7r° according to the distribution measured in Ref. [231] • All other 
1]' decays are generated using phase space. For each decay mode, we try all combinations 
for assignment of the selected tracks and showers (although in practice more than a single 
successful one per event is rare). 

We also tally an event satisfying the selections if instead of lying inside the light meson 
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signal mass window it has mass within one of the two sideband regions indicated in Table |TT1 
Similarly, the number of events with 7r"'"7r~ recoil mass inside the sidebands is accumulated 
for each J/ip mode. To allow for possible double-counting of non-J/ip and light meson 
mass-sideband backgrounds, we only subtract half of the window-size-scaled dipion recoil 
mass sideband number for each mode, and assign a systematic uncertainty of the full scaled 
number. 

We subtract feedacross attributable to any radiative J/ip or ip{2S) decays into 77?°, 777, 
and 777', including the decay modes not selected. To do so, we employ MC samples that 
are normalized to the observed net yields and r] ' relative branching fractions found here, 
except for i], cu, and 7r° decay modes not probed, for which PDG [l2| branching fractions 
are used. Continuum data counts satisfying the selections are subtracted after scaling by 
relative luminosity and 1/s. The mass sideband subtraction is performed in the MC and 
continuum samples, assuring that the efficiency, feedacross, and continuum backgrounds are 
computed correctly. This procedure implicitly assumes that backgrounds other than from 
continuum and feedacross are linear in pseudoscalar candidate mass; it explicitly accounts 
for nonlinear backgrounds from continuum and feedacross. 

Applying the above-described selection criteria to our data and MC samples yields the 
results shown in Table [U in which raw yields for signal and sidebands are given along with 
appropriately scaled feedacross and continuum corrections. The signal efficiency is also 
given for each mode, which ranges from 5% to 36%. All modes except ip{2S) 'yn^ and 
'?/'(3770) — i> 77r° have very small backgrounds. No statistically significant signals are seen for 
^(2^) ^ -ff], i:{2S) 77r0, or ^(3770) -> 771°. 

The light meson mass distributions for all modes are shown in Figs. [T}|31 plots for J/ip 
(other) modes have logarithmic (linear) vertical scales. Background from e~^e~ 37 that 
contaminates ^'(25') — > 77r° and ?/'(3770) — > 77r°, visible in Figs. [U^b) and|l](a), respectively, 
is irregular in shape but modeled well by the scaled continuum data. There is also a nonlinear 
background shape for modes involving rj 77r+7r^ decays, but it is well-modeled by the 
MC as it is due to feedacross from rj — *• 71~^71~tc^. A typical '?/'(25')-to-J/?/'-transition dipion 



recoil mass distribution appears in Ref. [22|, showing a tiny and flat non- J /ijj background; 
this is typical of all modes. Ample signals for J/ip —>■ 77r°, J/ip ^ jr], J/ip —>■ 777', and 
'4^{'2S) —>■ 777' are apparent and well- modeled by the MC-predicted shapes, aside from very 
small backgrounds that are approximately linear in mass. 

We make the flrst observations of e+e~ — > 777^''' at y/s = 3.773 GeV with statisti- 
cal signiflcances of 5.0(6.4)(T and observed cross sections (0.17l[5:[ji ± 0.03) pb (77/) and 
(0.2llo!o5 i 0.03) pb (77/'). The fraction of these values attributable to continuum cannot 
be estimated accurately from our ^/s = 3.671 GeV dataset due to its relatively small size: 
zero observed 77/^'^ events yields upper limits at 90% C.L. of < 2 pb for both proceses. How- 
ever, continuum cross sections for e~^e~ — 'yrj^'^ at y/s = 10.58 GeV have been measured by 



BaBar ^] to be 4.5li;i ± 0.3 (5.4 ± 0.8 ± 0.3) fb. Extrapolating [15|] to the charmonium 
energy region requires knowledge of the dependence of the form factors jHI F (,) (g^) for 



7* — s> 777*^'^ According to CLEO measurements [iJ] of the 77* — rj^'^ spacelike form fac- 
tors, \q'^F^(,){q'^)\ is 15 ±10 (5 ± 5)% smaller than that at = 7?i^[T(4S')] for g^ in the 
charmonium region. Continuum cross sections for e^e~ — * 777^'^ are expected [2^ to scale 
as |F^w(g2)|2, leading to an estimate of 0.19 ± 0.07 (0.25 ± 0.05) pb at ^/s = 3.773 GeV. 
These levels indicate a dominant continuum component in our signals. Other potential 
sources of background for the '?/'(3770) data sample are decays from the tail of the tp{2S), 
e~^e~ — > ip{2S) 7P at -^i = 3.773 GeV, and radiative returns (rr) to the peak of the 
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il}{2S), e~^e~ 7rrV'(25') 7rr7-P- Using the methods of Ref. [26|, these contributions are 
estimated to be neghgible: the former due to its inherently small rate and the latter because 
such events will fail the kinematic selections due to the presence of 

Notable in Fig. [2](c), [2](e), and[2](g) is the absence of any ip{2S) 777 signal; we set an 
upper limit on the cross section at y/s = 3.686 GeV of < 1.1 pb, about one tenth of that 
expected if R2 = Ri were true. We can rule out destructive interference with a continuum 
7?7 signal as the primary cause of the deficit because, as discussed above, the continuum 
cross section is too small. 

The statistical errors on efficiencies and subtractions for sideband, feedacross, and con- 
tinuum are combined with the statistical uncertainties on event counts in the data. For 
systematic uncertainties associated with finding tracks and showers, we assign 0.3% per 
track and 0.4% per photon y/7,] on a mode- by-mode basis, accounting for correlations in 
such uncertainties in arithmetic computations involving two or more modes. For the looser 
photon criteria applied to '?/'(3770) modes, a larger uncertainty of 3% per photon is assigned. 
Systematic uncertainties of 1% (relative) are assigned to efficiencies, uncorrelated mode-to- 
mode, and account for any dataset or trigger modeling dependence. We use rj branching 
fractions and uncertainties from the PDG fit [l2|- As a systematic cross-check, we compare 
ratios of J/ip ^ 'yrj corrected yields in different decay modes and find consistency with our 
previous measurement [2l| and Ref. fl2| . 

The net event yield for each channel listed in Table [T] is divided by its respective effi- 
ciency, intermediate decay branching fraction (from PDG [l2| for r] and this analysis for 
7]'), and total number of charmonia produced. The resulting branching fractions appear 
in Table IllH where multiple r]^'^ sub-mode values have been combined in a weighted aver- 
age. For J/ip ^ 77r°, 777, 777', and 4'{2S) — >• ■jrj' the results are consistent with previous 
measurements. We set improved upper limits for ip{2S) — *• 771° and 777 and first upper 
limits for iIj{3770) — > 77?°, 777, and 77/'. These upper limits are computed using simulated 
trials in which Poisson pseudo-random numbers are thrown for background and signal lev- 
els, accounting for sideband-window-size scaling, systematic errors, and fiuctuations in the 



observed backgrounds, in a manner similar to that in Ref. [27 . 

The selected J/ip jr]' events allow measurement of branching fractions for the five 
most common rj ' decay modes. Ratios of branching fractions follow from the J /ip ^ "^rj' 
entries in Tableland the world average value for B{ri —>■ 77) fl3\. The resulting ratios 
appear in Table HVl and represent the most precise individual measurements [l2t]. Absolute 
branching fractions for these five modes can also be obtained by constraining their sum. 
After combining branching fractions for 7r+7r~7r° and Ti~^Tr~e~^e~ [281] with those for 37r° and 
fj'~^fi~'y [l3], (99.2 ± 0.2)% of all 7/' decays remain for the largest five decay modes. We 
impose this constraint, accounting for correlated systematic uncertainties, and obtain the 
absolute branching fractions in Table IIV[ The values shown are consistent with those from 
the PDG constrained fit 12(] and of comparable precision. 

In summary, we have described improved measurements of and new limits on branching 
fractions for charmonium decays to a photon and a light pseudoscalar meson. We also have 
performed the first simultaneous measurement of the five largest r]' branching fractions, 
attaining improved precisions. Our result for B{J/i/j — *• jTr^) improves pre cision and confirms 



earlier central values 12|, yielding a better experimental constraint [l5| on the 7* — > 77r° 
vertex in the timelike regime at = m^lJ/ip). For the ratio of 7/ : 77 ' production rates from 
each resonance we obtain Ri = (21.1 ±0.9)% for J/ip and R2 < 1.8% at 90% C.L. for V(25), 
where for both results statistical and systematic uncertainties from the input branching 
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fractions from Table III have been combined in quadrature after accounting for correlations. 
Such a small value of R2 is unanticipated, posing a challenge to our understanding of cc 
bound state decays. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. D. Cronin-Hennessy and A. Ryd thank the A. P. Sloan 
Foundation. This work was supported by the National Science Foundation, the U.S. De- 
partment of Energy, the Natural Sciences and Engineering Research Council of Canada, and 
the U.K. Science and Technology Facihties Council. 
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TABLE I: For each decay mode, the number of events in the signal region (N); the total number in 
sideband intervals, both for meson candidate mass (Sm) and, where applicable, transition dipion 
recoil mass (Sj), without scale factors applied; the sum of scaled continuum and MC feedacross 
backgrounds (F), which is negative when scaled sidebands exceed signal counts; and the MC 
efficiency (e). 



Mode 


iV 


Sm 


Sj 


F 


6(%) 


J/tp ^ 77r° 


113 


9 


33 


-2 


31.1 


77/(77) 


1137 


80 


69 


4 


26.3 


^ 7r/(37r0) 


256 


21 


6 





7.0 


— ^ 7?7(7r"'"7r~7r^) 


217 


5 


30 





9.0 


— > 7r/(77r"'"7r~) 


105 


62 


12 


-12 


21.0 


7r/'[7r+7r~7?(77)] 


1208 


36 


56 


3 


14.2 


77/ '[vr'^TT*' 77(77)] 


245 


63 


16 


-2 


4.9 


— > 777'(77r+7r^) 


3205 


133 


129 


5 


21.8 


— ^ 77/'[7a;(7r'^7r~7r°)] 


71 


18 


8 





5.6 


111' ill) 


317 


230 


16 





22.7 


7^(25) ^ 77r° 


31 


89 


_ 


-33 


20.0 




2 


4 


_ 





15.8 


— ^ 777(7r"'"7r~7r^) 


1 


1 


_ 





12.6 


— 77/(77r+7r~) 


1 


7 




-1 


32.1 


— > 77/'[7r"'"7r~7?(77)] 


120 


3 







22.6 


77/ '['/r'^Tr'^ 7/(77)] 


46 


39 







11.2 


^ 77/'(77r"'"7r~) 


343 


91 




1 


35.3 


— > 77/'[7a;(7r"'"7r~7r°)] 


12 


6 




-1 


9.9 


^^(3770) ^ 77r° 


331 


1396 




-468 


19.8 


^ 7r?(37r0) 


9 










19.4 


^ 77/(7r"''7r~7r'^) 


7 


1 







23.2 


— > 777'[7r"'"7r~77(77)] 


8 


1 







33.9 


— ^ 77/'[7r"'"7r~7;(37r°)] 


3 










7.3 


— ^ 77/'[7r"'"7r~r/(7r+7r~7r'^)] 


5 










13.5 
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TABLE II: Mass windows and sideband (SB) intervals in MeV. 



Mode Window Lower SB Upper SB 



7r+ 


TT -recoil 


3087-3107 


2980-3080 


3114-3214 




77 


110-160 


50-100 


170-220 


77 - 


-> 77r^7r~ 


535-560 


460-510 


585-635 




— > other modes 


500-580 


400-480 


600-680 


77' 




920-995 


730-880 


1035-1185 




other modes 


945-970 


890-940 


975-1025 


UJ - 


7r+7r~7r'' 


750-814 







TABLE III: Branching fractions from this analysis, showing statistical and systematic uncertainties, 



respectively, and PDG I4]. The rightmost column shows the difference between the two in units 
of standard error (a). Upper limits are quoted at 90% C.L. Entries in the last two rows include the 
effects of estimated continuum background and ignore (include) maximal destructive interference 
between '(/;(3770) and continuum sources. 



Mode This result (10"^) PDG (10"^) #(T 

J/V' ^ 7vrO 0.363 ±0.036 ±0.013 0.33+^;^^ oT 

^77/ 11.01 ±0.29 ±0.22 9.8 ±1.0 1.2 

-^jrj' 52.4 ±1.2 ±1.1 47.1 ±2.7 1.7 

ip{2S) 77r° < 0.05 < 54 

^77/ < 0.02 < 0.9 

-^'jrj' 1.19 ±0.08 ±0.03 1.36 ±0.24 -0.7 

^^(3770) ^ 77r° < 2 - - 

-^'jrj < 0.2 (1.5) 

-^'jrj' < 0.2 (1.8) 



TABLE IV: Relative and absolute branching fractions for 77' decays from this analysis, showing 
statistical and systematic uncertainties. 



Mode B/Bjir+TT-T]) (%) B{%) 

vr+vr-r/ = 100 42.4 ±1.1 ±0.4 

7r°7r°r/ 55.5 ± 4.3 ± 1.3 23.5 ± 1.3 ± 0.4 

■y-K+TT- 67.7 ±2.4 ±1.1 28.7 ±0.7 ±0.4 

70; 5.5 ±0.7 ±0.1 2.34 ±0.30 ±0.04 

77 5.3 ±0.4 ±0.1 2.25 ±0.16 ±0.03 



10 



2540309-004 




M (yy) (GeV) 

FIG. 1: Two-photon mass distributions for jir^ final states from (a) J/ip and (b) ip{2S). Points 
show the ip{2S) data; dashed histogram, the luminosity-scaled continuum data; solid line histogram, 
signal MC; and dotted histogram, MC feedacross from other 7P decays. Upper solid arrows show 
nominal selection criteria; all selection criteria are in place except those upon plotted mass. Lower 
dashed arrows show sidebands. 
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FIG. 2: Distributions of rj candidate mass for 777 final states from J/'i/; (left) and ip{2S) (right). 
The r] decay modes are (a) 77; (b), (c) 3tt^; (d), (e) vr+Tr^vr'^; and (f), (g), 77r'^7r~. Symbols are 
defined in Fig. 1. 
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FIG. 3: Distributions of rj' candidate mass for jrj' final states from J/ip (left) and ip{2S) 
(right). The 7]' decay modes are (a), (b) 7r+7r~Ty(77); (c), (d) tt^tt'' 7/(77); (s), (f) 7vr+7r~; (g), 
(h) 7u;(7r^7r~7r'^); and (i) 77. Symbols are defined in Fig. 1. In (c), (d), and (i), only parts of the 
sideband intervals are shown. 
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FIG. 4: Distributions of light quark meson candidate mass for y/s = 3.773 GeV final states of (a) 
77r'^; (c), (e) jr]; and (b), (d), (f) 717'. Decay modes are (b) 7?7'[7r+7r^77(77)], (c) 7?7(37r'^), (d) 
7?7'[7r'^7r~r/(3-7r'^)], (e) 7r/(7r^7r~7r'^), and (f) 7?7'[7r"^7r~?7(7r^7r~7r'^)]. Symbols are defined in Fig. 1. 
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